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ABSTRACT 



The origin of color gradients in elliptical galaxies is examined by comparing 
model gradients with those observed in the Hubble Deep Field. The models are 
constructed so as to reproduce color gradients in local elliptical galaxies either 
by a metallicity gradient or by an age gradient. By looking back a sequence 
of color gradient as a function of redshift, the age - metallicity degeneracy is 
solved. The observed color gradients in elliptical galaxies at2; = 0.1 — 1.0 agree 
excellently with those predicted by the metallicity gradient, while they deviates 
significantly from those predicted by the age gradient even at z ^ 0.3 and the 
deviation is getting larger with increasing redshift. This result does not depend 
on cosmological parameters and parameters for an evolutionary model of galaxy 
within a reasonable range. Thus our results clearly indicate that the origin of 
color gradients is not the age but the stellar metallicity. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: 
evolution — galaxies: formation 
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1. INTRODUCTION 

Stellar populations in an elliptical galaxy are not uniform. Stars at a galaxy center are 
redder than those in the outer region and colors in a galaxy are progressively getting bluer 
with an increasing radius (e.g., Vader et al. 1988; Franx, lUingworth, & Heckman 1989; 
Peletier et al. 1990a; Peletier, Valentijn, & Jameson 1990b). Peletier et al. (1990a) made 
surface photometry in the U-, B-, and i?-bands for a sample of 39 nearby elliptical galaxies 
and examined their color gradients in U — R and B — R. They found that typical color 
gradients A(C/ — i?)/Alogr and A(i3 — i?)/Alogr are —0.20 mag/dex and —0.09 mag/dex, 
respectively, and demonstrated that a dispersion of the color gradients is small, i.e., only 
0.02 mag /dex in both colors. Since many of elliptical galaxies show radial gradients in 
line strengths such as Mg2, Fei(5270 A) and Fe2(5335 A) (e.g., CaroUo, Danziger, & Buson 
1993; Davies, Sadler, & Peletier 1993; Gonzalez 1993; Kobayashi & Arimoto 1999), the 
color gradients have been naively assumed to originate from a metallicity gradient inside a 
galaxy. 

However, such an interpretation for the origin of the color gradient is premature, 
because stellar populations of either higher metallicity or older age can make a galaxy 
redder. This problem, which is called the age-metallicity degeneracy, was first pointed out 
by Worthey, Trager, & Faber (1996) and then discussed by Arimoto (1996). For example, 
the degeneracy makes it difficult to interpret the origin of the tight correlation between 
colors and magnitudes of elliptical galaxies; brighter elliptical galaxies tend to have redder 
colors. This correlation called color-magnitude (CM) relation can be excellently reproduced 
by a metallicity sequence with a galactic wind model based on a monolithic collapse scenario 
(e.g., Arimoto & Yoshii 1987), where more massive elliptical galaxies should be enriched 
more in metals and thus become redder. However, Worthey et al. (1996) claimed that an 
age sequence of eUiptical galaxies can equivalently reproduce the CM relation if brighter 
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elliptical galaxies are older and thus redder. To break this degeneracy, Kodama & Arimoto 
(1997) built up two model sequences (metallicity and age sequences) which are normalized 
to reproduce the CM relation of ellipticals in Coma cluster by using their evolutionary 
synthesis model. They compared the evolution of the model CM relation with the observed 
relations of eUipticals in distant clusters. The CM relations produced by a metaUicity 
sequence agree with the observed relations out to z ^ 1, while those produced by an age 
sequence deviates from the observed ones significantly even at 2; ~ 0.2 — 0.3, showing the 
origin of the CM relation is primarily a metaUicity variation with a galaxy mass. The CM 
relation can also be reproduced in a hierarchical galaxy formation scenario (Kauffmann & 
Chariot 1998). Although a model of Kauffmann & Chariot (1998) allows a more extended 
period of star formation in elliptical galaxies, it shows that the CM relation is produced by 
a metallicity variation. It is thus worth emphasizing that the interpretation of CM relation 
with a metallicity sequence is robust and independent of detailed assumptions on galaxy 
formation processes. 

Recently, the origin of color and line strength gradients in elliptical galaxies has been 
discussed with elaborate models. Martinelli, Matteucci, & Colafrancesco (1998) tried to 
reproduce the color gradients by assuming that a galactic wind blows later in the inner 
part of a galaxy due to a deeper potential well defined mainly by dark matter. Adopting a 
multi-zone model that takes into account gas dynamics, local star formation, and chemical 
evolution, Tantalo et al. (1998) reproduced radial gradients of colors and hue strengths to 
some extent. Nevertheless, these detailed modelings are not fully successful; the inner part 
of a galaxy becomes too iron enriched due to an extended period of star formation there. 

In this paper, adopting a much simpler approach without entering into details of 
physical processes of galaxy formation and evolution, we try to depict essential aspect of the 
origin of color gradients. Our approach is similar to that adopted by Kodama & Arimoto 
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(1997) for studying the origin of the CM relation. By using a population synthesis model, 
we first make two different model galaxies, both of which can reproduce a typical color 
gradient of elliptical galaxies ai z = 0, by changing either mean stellar metallicity or age, 
and let evolve them back in time. The evolution of color gradients thus predicted are then 
confronted with the observed ones in distant ellipticals extracted from the Hubble Deep 
Field North (HDF-N; Williams et al. 1996). As was demonstrated by Kodama & Arimoto 
(1997), the best way to disentangle the age and metallicity effects on galaxy colors is to 
look back galaxies at high redshift. 

It should be noted that dust extinction in elliptical galaxies may have some effects 
on the color gradients (Goudfrooij & de Jong 1994; Wise & Silva 1996). Witt, Thronson, 
& Capuano (1992) calculated a radiative transfer within elliptical galaxies by assuming 
a diffuse distribution of dust and suggested that surface brightness profiles and color 
gradients could be well reproduced by dust effects. In fact, it has been beheved that the 
IRAS detected FIR emission for about half of the elliptical galaxies observed. However, 
since many of the detections were made with a ~ 3o" threshold, a significant fraction of the 
previously claimed detections may be spurious. Only 12 — 17% of the observed elliptical 
galaxies are detected with a sufficiently high confidence level (Bregman et al. 1998). 
Therefore, in this paper, we have chosen to focus on age and metallicity effects only. Effects 
of dust extinction on the color gradients are still open and will be studied in our subsequent 
paper. 

This paper is organized as follows. The model prescription is given in § 2. The sample 
selection and data reduction of elliptical galaxies in the HDF-N are described in § 3. The 
observed color gradients are compared with the models in § 4. Discussion and conclusions 
are given in § 5 and § 6, respectively. The cosmological parameters adopted throughout 
this paper are Hq — 50 km s~^ Mpc~^, flo — 0.2 and A = unless otherwise noted. 
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2. MODELS 

The color gradients can be reproduced by either a metaUicity gradient or an age 
gradient, but these effects are degenerate ai z = 0. We attempt to break up such degeneracy 
by confronting model color gradients with the observed ones of ellipticals at high redshifts 
in the HDF-N. We have built two sequences of evolutionary models under the alternative 
assumptions: 1) the color gradients originate from the metaUicity gradient of old stellar 
populations, or 2) the color gradients arise from the age gradient of stars which have the 
same metaUicity within the galaxy. 

The metaUicity sequence is constructed as follows: We assume that the star formation 
lasted progressively longer towards galaxy center; i.e., the galactic wind blew later in the 
inner region, so that the mean stellar metaUicity becomes higher and the colors become 
redder. Even at the galaxy center, however, the star formation stopped at early times 
(tgw = 0.83 Gyr) and the stellar ages are almost the same at everywhere within a galaxy; ~ 
15 Gyr at 2; = (see Tables 1 and 2). The assumption seems to be justified, because the 
duration of star forming activity is determined from the local dynamical potential (Larson 
1974) and should be longer in the inner region. To build a metaUicity gradient, however, 
this hypothesis is not unique; other hypotheses such as higher star formation efiiciency, a 
larger yield (i.e., a fiatter initial mass function (IMF) slope), or higher metaUicity of infalling 
gas towards the galaxy center can also produce the metaUicity gradient. However, such 
details are not essential in the present study, because so far as the mean stellar metaUicity 
increases in the inner part of a galaxy while stars are uniformly old, the resulting evolution 
of color gradient is essentially the same. We should stress this, because this paper aims for 
showing how color gradients of stellar populations behave as a function of lookback time, 
not for seeking the best model that can explain the observed color gradients of sampled 
galaxies in the HDF-N. The galactic wind is conventionally introduced to build up the 
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metallicity gradient and it might be possible that such gradients can also be established if 
elliptical galaxies formed via hierarchical clustering. 

The age sequence is constructed as follows: The star formation started earher in the 
inner regions, but lasted equally long (tg^ = 0.83 Gyr) everywhere in a galaxy. The period 
chosen here is the same as that assumed for the center of metallicity sequence model. In 
this way, the mean stellar ages t^ge become older and colors redder in the inner regions, 
while the mean stellar metaUicities are almost the same everywhere (see Tables 1 and 2). 
It is hard to comprehend the physical process to produce such age gradient, although it 
may occur if a series of young dwarf galaxies accrete onto a massive galaxy and are tidally 
disrupted at the outer region of a galaxy before falling into the center. However, it is not 
our aim anyway to construct a physically motivated age sequence model. 

Although we do not discuss in this paper in detail, we have tried to build an alternative 
age sequence model in which star formation began everjrwhere 15 Gyr ago and lasted longer 
at the outer parts of the galaxy, keeping the same metallicity everywhere as that at the 
galaxy center. As a result, the mean stellar ages become progressively younger from the 
galaxy center towards the outer region. However, this model fails to reproduce the observed 
color gradients in nearby ellipticals. If old stellar populations arc contained at the outer 
regions and if the mean metallicity is as rich as that at the center, there is no way by 
making their mean ages young to build up the color gradients as steep as the observed ones. 

Both metallicity and age sequences are constructed to reproduce typical color gradients 
observed for nearby elliptical galaxies; A{B — i?)/Alogr = —0.09 ± 0.02 mag/dex (i.e., 
—0.07 and —0.11 mag/dex) and B — R = 1.633 mag at r = re/10 (re refers to an effective 
radius) which is derived from the mean B — R color at re/2 and A(i^ — i?)/Alogr = —0.09 
mag/dex (Peletier et al. 1990a). We also construct the models which reproduce the color 
gradients of —0.09 ± 0.04 mag/dex (i.e., —0.05 and —0.13 mag/dex) within which the 



-8- 



gradients of almost all the sample ellipticals in Peletier et al. (1990a) are included. Thus we 
have a set of four model sequences each for the age sequence and the metallicity sequence. 

To calculate gradients of photo-chemical properties of a galaxy, we construct a 
galaxy model consisting of a 'shell' at each radius, and assume that each 'shell' evolves 
independently. The star formation rate (SFR) in a 'shell' is proportional to the gas fraction 
(Schmidt law) with a time scale of 0.1 Gyr and a primordial gas is supplied to the shell 
with a rate of exp(— i/0.1 Gyr). An infall of a primordial gas onto the shell may sound 
rather strange, since it is more likely that the inner shells suffer from the infall of enriched 
gas from the outer shells. However, a proper modeling of such effect is beyond our scope. 
Since the infall rate we employed in the present study is significantly high, the resulting 
evolutionary behavior of the color gradient remains almost the same even if we adopt the 
'simple' model prescription instead of the infall model. For an IMF, a power-law mass 
spectrum with a slope oi x — 1.10 is assumed in the range of O.OSM© < M < SOM©. We 
note that this IMF slope is the same as in Kodama et al. (1998), and was introduced 
to increase chemical yields to reproduce the reddest end of the Coma CM relation. The 
nucleosynthesis yields of SNe la and SNe II are taken from Tsujimoto et al. (1995). The 
chemical enrichment by SNe la is calculated with a metaUicity dependent SN la rate, of 
which detail formulation is given by Kobayashi, Tsujimoto & Nomoto (1999). The spectral 
evolution is calculated using a spectral synthesis database of Kodama & Arimoto (1997). 
The adopted parameter values for galactic wind epoch and stellar age t^ge at each radius 
for the four color gradients are hsted in Table 1. The resulting photo-chemical properties 
of a galaxy aX z = are summarized in Table 2, including gas abundances and colors at 
each radius. It should be noted that our models have an SN Il-like abundance pattern, 
though the value [Mg/Fe] ~ -1-0.4 is slightly larger than that suggested by the observational 
estimates (Worthey, Faber & Gonzalez 1992; Kobayashi & Arimoto 1999); a larger [Mg/Fe] 
ratio is a result of the assumed flat IMF {x = 1.10). 
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3. GALAXIES FOR COMPARISON 

3.1. Sample Elliptical Galaxies 

To compare theoretical color gradients with those of elliptical galaxies at high redshifts, 
the archival data of the HDF-N (WiUiams et al. 1996) are used. Our sample galaxies are 
selected from those brighter than Isu,ab = 22 mag, in such a way that we can derive 
reliable surface brightness profiles and color gradients. All these galaxies have spectroscopic 
redshifts. EUiptical galaxies are identified by using a bulge-to-total (B/T) luminosity ratio 
derived by Marleau & Simard (1998) who obtained the B/T ratios in the /814-band of HDF 
galaxies brighter than /814 = 26 mag by decomposing quantitatively the surface brightness 
profile into the bulge and the disk components. We define galaxies with B/T > 0.5 as 
ellipticals according to Marleau & Simard (1998). Our resulting sample consists of ten 
elliptical galaxies with redshifts spanning from z = 0.089 to 1.015 as listed in Table 3. It 
is important to note here that our sample selection does not rely on any color information. 
Franceschini et al. (1998) also made a sample of elliptical galaxies in the HDF using surface 
brightness profiles. Their selection, however, is not based on the bulge - disk decomposition. 
Consequently, within /814,ab < 22 mag, our sample is smaller than that of Franceschini 
et al. (1998), except for the galaxy 4-241.1 a,t z — 0.318. The Franceschini et al's sample 
includes several galaxies having the B/T ratio smaller than 0.5, which are presumably disk 
dominant galaxies and thus are not included in our sample (see also Kodama, Bower & 
Bell 1999). It should be noted that B/T > 0.5 in the /su-band may be slightly loose to 
isolate eUiptical galaxies. B/T > 0.5 in the S-band roughly corresponds to local early- type 
galaxies (e.g., Simien & de Vaucouleurs 1986). By using bulge and disk models of Kodama 
et al. (1999), we found that the /814-band B/T ratios of ellipticals are larger than ~ 0.7 
at 2; = and change little from 2; = to 1 in the observer's frame. Thus our sample may 
include galaxies as late as Sa, though they must be minority because most of the B/T 
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ratios of our sample galaxies are larger than 0.7 (Table 3). 

We measured Veoe — -^8i4 colors of our sample galaxies and present them in Table 3 and 
Figure 1. These colors are obtained for a 10 kpc aperture. Seven of the sample galaxies 
have red colors consistent with those for eUipticals with passive evolution as shown in 
Figure 1 (hereafter we call these galaxies 'red ellipticals'), and the other three galaxies 
(4-241.1, 2-251.0, and 4-928.0) have blue colors (see also Figure 2; hereafter referred to 
as 'blue eUipticals' for convenience). Azimuthally averaged surface brightness profiles of 
the sample ellipticals in the /814-band are shown in Figure 2. The average is taken along 
an ellipse fitted to an isophote of the /814-band image. These profiles are well represented 
by an r^/^ law for both the red eUipticals and the blue eUipticals. Effective radii (re) of 
the sample eUipticals are obtained by the r^/^ fitting and results are shown in Table 3. 
The fitting is done by removing data points in the inner and outer regions, since they are 
unreliable due to an effect of a point spread function (PSF) and low signal-to-noise ratios 
(S/Ns), respectively. The obtained effective radii range from Ve —lA to 7.9 kpc, which are 
typical for nearby giant ellipticals (e.g.. Bender, Burstein, & Faber 1992). 



3.2. Color Gradients 

To examine color gradients of our sample galaxies, the images in F606W (Veoe)- and 
F814W (/8i4)-bands were used, because S/Ns of the images are better in these bands than 
in the other bands (F300W and F450W). First, the sky value around each object was 
determined by using the "phot" task in the IRAF apphot package ("mode" in an annulus 
with an inner radius of mostly 3" ~ 4" and a width of 0."4 is adopted) and was subtracted. 
Next, the angular resolutions of the blue images and the red images were adjusted and a 
Veoe — hu color map of each galaxy was made. Figure 3 shows the resulting color maps. 
Most of the color maps for the red ellipticals show symmetric structure and ordinary color 
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gradient. In two cases (2-456.0 and 2-121.0), slight asymmetry in the color distribution 
is seen. This does not seem to be caused by a misalignment of the images. For the blue 
ellipticals, variety of the color distribution is seen, which we will discuss in §5.1. Next, a 
1^06 — -^814 color profile is derived with the azimuthally averaged radial profiles as a function 
of semi-major radius from the galaxy center. Resulting color profiles are shown in Figures 
4 and 5 for the red ellipticals and in Figure 7 for the blue ellipticals. It should be noted 
that in plotting observed color profiles, a zero-point of the observed color for each galaxy 
is shifted to be compared with those by the models. The amounts of these shifts are given 
in the last column of Table 3. Most of the shifts are less than 0.2 mag; for intrinsically 
luminous objects the colors tend to be shifted to bluer colors and vice versa. Thus the color 
differences are presumably caused by an effect of the color-magnitude relation of elliptical 
galaxies. The zero-point offsets could also come from photometric errors in HST data 
(e.g., Holtzman et al. 1995; Ellis et al. 1997; Kodama et al. 1998) and from uncertainty 
of the population synthesis model. Needless to say, zero-point offsets do not affect the 
color gradients. An error bar attached to each data point in the color profile includes 
an photometric error, a local sky subtraction error, and a dispersion of colors along each 
elliptical isophote. The data points of 2-456.0 and 2-121.0 have relatively large error bars in 
the middle region of the profiles; these objects have asymmetric color distributions which 
results in the relatively large dispersions of colors in the regions. 

Finally, slopes of the color profiles A(V606 — -^8i4)/'Alog(r/re) of the sample galaxies are 
derived by applying a least square fit to the color profiles. This fit was done after removing 
the unreliable data points in the outer region (r > re) and in the innermost region. The 
removed data points in the inner most regions are shown by crosses in the color profiles 
(Figures 4 and 5); they clearly deviate from the r^l^ fits presumably due to the effect of a 
PSF as seen in Figure 2. The slopes of the color profiles for the blue eUipticals presented in 
Figure 7 are not derived. 



4. ORIGIN OF COLOR GRADIENTS 



Once the spectra are calculated at various radii of a model galaxy in such a way that 
they reproduce the color gradient A(B — R)/ /S.\og{r/r^) = —0.09 mag/dex at 2; = 0, the 
model gradient in any colors at any redshifts can be predicted to compare with the observed 
gradients. The model color gradients at the redshift of each sample eUiptical galaxy are 
over-plotted in each panel of Figures 4 and 5 for the age gradient and the metaUicity 
gradient models, respectively. Solid and dotted lines in each panel show the gradients in 
Veoe — -^814 color in the observer's frame corresponding to the gradients of —0.09 ± 0.02 
mag/dex and —0.09 ± 0.04 mag/dex in the B — R color at 2; = 0, respectively. 

It is clearly shown that the model gradient made by the age gradient begins to deviate 
from the observed ones at a redshift of 2; ~ 0.3 and the deviation is getting worse as a 
redshift increases (Figure 4). On the contrary, the model gradient made by the metaUicity 
gradient agrees well with the observed color gradients within the effective radius at any 
redshift from 2; ~ 0.1 to ~ 1 (Figure 5). Considering that the model gradients are calibrated 
only at 2; = 0, we insist that the agreement between the model and the observed gradients 
is excellent. Such agreement is more clearly seen in Figure 6 which shows both the observed 
color gradients and the model gradients as a function of redshift. Solid and dotted lines 
show the evolution of model color gradients in Veoe — in the observer's frame for the 
metaUicity gradient and the age gradient, respectively. These model gradients are calculated 
in a region of — 1 < log r/re < for the metaUicity gradient and for the age gradient at 
z < 0.5, and —1 < logr/re < —0.5 for the age gradient at 2; > 0.5. Dots show the observed 
gradients of the sample galaxies within calculated by a least square fitting to the data 
points shown as filled circles in Figures 4 and 5. Error bars indicate the fitting errors of la. 



5. DISCUSSION 
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5.1. Blue Galaxies in Our Sample 

Since our sample is selected independently of galaxy colors, relatively blue elliptical-like 
galaxies are included in our sample (3/10). Among these three, 4-241.1 (z = 0.318) has 
a nearly flat color proflle with a slightly bluer color in the inner part (Figure 3 and the 
top panel of Figure 7). The spectrum of this galaxy (from the Hawaii Active Catalog; 
Cowie 1997) shows Balmer absorption hues shortward of 4000A in the rest frame as well 
as emission lines such as [OII]A3727, B.(3, [OIII]A4959, A 5007, and Ha. Thus this galaxy 
is expected to have some young stellar populations. The galaxy 2-251.0 {z — 0.960) has 
a color gradient in the opposite sense to those of ordinary ellipticals; i.e., the inner region 
of the galaxy is bluer than the outer region (Figure 3 and the middle panel of Figure 
7). Since this galaxy has been suggested as an AGN (Franceschini et al. 1998), this 
opposite gradient is probably caused by AGN and/or associating star burst activity in the 
central part of the galaxy. Although both galaxies seem to have the younger ages, our 
model by the age gradient cannot reproduce their color gradients. Accordingly, neither the 
metallicity gradient nor the age gradient can explain their observed color gradients. These 
galaxies presumably have on-going star formation (or post-starburst) in the inner part of 
the galaxies, which may be caused by a galaxy interaction as expected from their rather 
disturbed morphology. Our model cannot apply to these cases. The other galaxy 4-928.0 
{z — 1.015) has a rather steep color gradient as seen in Figure 3 and the bottom panel of 
Figure 7. The observed color profile of the galaxy is shifted by 0.4 mag to be compared 
with the model color gradients (solid and dotted lines). The observed color gradient of this 
galaxy may be explained by the age gradient if we tune ages of stellar populations in the 
model galaxy. However, the mean age is set to be 0.2 Gyr at re and 9 Gyr at the innermost 
point. Such a large age at the center cannot be realized at this redshift {z = 1.015) unless 
a different set of cosmological parameters is adopted. Finally, note that colors of these 
galaxies are not consistent with those of passively evolving ellipticals formed at high z. 
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Much larger sample is needed in future works to find whether the population of the blue 
galaxies is really the minority in the whole population of elliptical galaxies. 

5.2. Gradients of Metal Absorption Line Strengths 

Kobayashi & Arimoto (1999) showed that central velocity dispersions of elliptical 
galaxies correlate well with their mean metallicities. They found, however, that this 
correlation has a significant scatter. They also showed that Mg2 gradients do not correlate 
with any physical parameters including a galaxy mass. Such a large dispersion and an 
absence of the correlation are contrary to what monolithic collapse simulations predicted. 
As the common origin for the large dispersion and the absence of the correlation, they 
suggested the following possibilities: (1) a significant age spread among elliptical galaxies, 
(2) age gradients inside elliptical galaxies, (3) dust obscuration, (4) effects of merging 
on metallicity gradients inside elliptical galaxies, and (5) intrinsic scatter which elliptical 
galaxies have had since they formed. Since we confirm in this paper that the color gradients 
are caused by the metallicity gradient, we can reject the possibihty (2). 

There seems to be a problem that color gradients in eUiptical galaxies do not correlate 
with their Mg2 gradients, though only a small number of elliptical galaxies are examined 
both for color gradient and Mg2 gradient (Peletier 1989). However, the absence of this 
correlation may be a spurious result caused by a different way to derive the gradients 
between the color and Mg2 gradients. Some effects of dust extinction and/or, possibly, 
observational errors on both gradients might suppress the true correlation. In any case, this 
problem will be investigated in detail in our forthcoming paper. 

Worthey et al. (1992) suggested that the color-magnitude relation of elliptical galaxies 
may originate from the fact that more luminous ellipticals tend to contain more Mg 



- 15 - 



relative to Fe. One can claim that it could be possible that the same holds for the internal 
gradients. However, Kobayashi & Arimoto (1999) have recently carried out a careful study 
of line strength gradients of 80 elliptical galaxies having the best quality of spectra but find 
no evidence supporting for the [Mg/Fe] gradient. Therefore, we conclude that the color 
gradients are not originating from the gradient of the [Mg/Fe] ratio. 

6. CONCLUSIONS 

The origin of color gradients in elliptical galaxies is examined. A typical color gradient 
of nearby elliptical galaxies is reproduced by two alternative model sequences, metallicity 
gradient and age gradient. This age-metaUicity degeneracy is broken by looking back the 
evolution of color gradient towards high redshifts; the predicted color gradients at high 
redshifts arc compared with the observed color gradients of ten elliptical galaxies in the 
HDF-N out to ^ ~ 1. We find that the observed color gradients of the seven red galaxies, 
whose colors are consistent with those of passively evolving galaxies, are in excellent 
agreement with the metallicity gradient at any redshift. These conclusions are independent 
of cosmological parameters and parameters for an evolutionary model of galaxy (IMF and 
SFR) within a reasonable range. Thus we conclude that the primary origin for the color 
gradient in elliptical galaxies is the metallicity gradient in the old stellar populations. 

This work was financially supported in part by Grant-in- Aids for the Scientific Research 
(No. 0940311 and 09740173) by the Japanese Ministry of Education, Culture, Sports and 
Science. CK and TK thank Research Fellowships of the Japan Society for the Promotion of 
Science (JSPS) for Young Scientists. 
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Table 1. Input Parameters 



log r/re 


-1 


-0.875 -0.75 


-0.625 


-0.5 


-0.375 


-0.25 


-0.125 











Metallicity Sequence {ta 


,ge = 15 Gyr) 










A(S-i?)/Alogr 








igw [Gyr] 










-0.05 


0.83 


0.74 0.64 


0.58 


0.53 


0.49 


0.45 


0.42 


0.39 


-0.07 


0.83 


0.69 0.59 


0.52 


0.47 


0.42 


0.38 


0.35 


0.32 


-0.11 


0.83 


0.63 0.52 


0.43 


0.37 


0.32 


0.29 


0.26 


0.23 


-0.13 


0.83 


0.60 0.48 


0.39 


0.33 


0.29 


0.26 


0.22 


0.20 






Age Sequence (tgw = 


0.83 Gyr) 










A{B-R)/Alogr 








W [Gyr] 










-0.05 


15.00 


14.35 13.65 


13.00 


12.40 


11.20 


10.20 


9.75 


9.20 


-0.07 


15.00 


14.00 13.10 


12.20 


10.60 


9.75 


9.00 


8.60 


7.80 


-0.11 


15.00 


13.60 12.00 


9.95 


8.90 


7.85 


7.18 


6.70 


6.20 


-0.13 


15.00 


13.35 11.10 


9.35 


8.35 


7.20 


6.65 


6.10 


5.63 



Note. — The second and third rows are the models constructed for the typical color gradients of 
A{B-R)/A logr = -0.09±0.02 mag/dex. The models in the first and fourth rows are for A{B-R)/A logr = 
—0.09 ± 0.04 mag/dex, within which gradients of almost all ellipticals in Peletier et al. (1990a) are included. 
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Table 2. Properties of Model Galaxies at 2; = 



log r/re 


-1 


-0.875 


-0.75 


-0.625 


-0.5 


-0.375 


-0.25 


-0.125 











Metallicity Sequence (fage 


— J^o ^yrj 










\og Z/Zq 


0.254 


U.Z4U 




n oni 
U.zUi 


0.182 


0.168 


0.142 


0.123 


0.101 




0.254 


U.zoU 


U.ZU4 


n 1 77 

U.i M 


0.153 


0.123 


0.092 


0.066 


0.034 




0.254 


n 01 A 
U.zlo 


U.i / / 


U.ioo 


0.084 


0.034 


-0.002 


-0.046 


-0.080 




0.254 


n 01 n 


U.iOc5 


U.iUi 


0.056 


-0.002 


-0.046 


-0.118 


-0.140 


[Fe/H] 


-0.069 


-U.Uoo 


n 1 1 


-U.ioi 


-0.151 


-0.165 


-0.192 


-0.213 


-0.236 




-0.069 


-U.iUU 




-U.ioo 


-0.181 


-0.213 


-0.245 


-0.273 


-0.307 




-0.069 


n 1 1 ^ 
-U. i io 


-U.ioo 


-U. iyy 


-0.254 


-0.307 


-0.346 


-0.394 


-0.431 




-0.069 


-U.izi 


n 1 71^ 
-U.i /O 


-U.ZOD 


-0.284 


-0.346 


-0.394 


-0.474 


-0.497 


[Mg/Fe] 


0.467 


U.4 / 


/I 71^ 
u.4 ( 


n /1 71^ 
U.4 / 


0.476 


0.477 


0.478 


0.479 


0.481 




0.467 


U.4/4 


n /1 71^ 
U.4 i 


n /1 7fi 


0.478 


0.479 


0.481 


0.483 


0.486 




0.467 


U.4 / 


n /1 7f^ 
U.4 ( 


n /1 70 
U.4 i y 


0.482 


0.486 


0.489 


0.493 


0.496 




0.467 


0.475 


0.477 


0.481 


0.484 


0.489 


0.493 


0.500 


0.503 


TT — P 
u — JX 




2.158 


2.144 


2.132 


9 118 


9 IDQ 


9 nsQ 




9 rififi 




2.166 


2.152 


2.134 


2.115 


2.098 


2.076 


2.055 


2.036 


2.014 




2.166 


2.143 


2.115 


2.085 


2.049 


2.014 


1.988 


1.957 


1.933 




2.166 


2.139 


2.101 


2.060 


2.029 


1.988 


1.957 


1.904 


1.888 


B-R 


1.632 


1.627 


1.620 


1.614 


1.608 


1.604 


1.595 


1.589 


1.583 




1.632 


1.624 


1.615 


1.606 


1.599 


1.589 


1.580 


1.572 


1.563 




1.632 


1.619 


1.606 


1.593 


1.578 


1.563 


1.552 


1.539 


1.528 




1.632 


1.617 


1.600 


1.583 


1.569 


1.552 


1.539 


1.516 


1.509 


Veoe — hi4 


0.576 


0.574 


0.572 


0.569 


0.567 


0.565 


0.562 


0.559 


0.556 




0.576 


0.573 


0.570 


0.566 


0.563 


0.559 


0.555 


0.552 


0.548 




0.576 


0.571 


0.566 


0.561 


0.554 


0.548 


0.543 


0.536 


0.531 




0.576 


0.570 


0.564 


0.556 


0.550 


0.543 


0.536 


0.525 


0.522 
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Table 2 — Continued 



loRr/;v 


-1 


-n.<S7r) 


-0.7o 


-0.()2."') 


-O.o 


-0.375 


-0.20 


-0.120 









Arc 




(7„„, = 0.8:-! G\-r) 










t'age 


10. uu 


14. OO 


lo.oo 


13.00 


12.40 


1 1 9n 

1 1 .ZU 


1 n 9n 

lU.ZU 


Q 7^ 

y. / 


Q 9n 
y.zu 




1 nn 


1 A nn 


1 1 n 


12.20 


10.60 


Q 7ii 
y. 1 o 


Q nn 
y .UU 


8 fin 

O.OU 


7 sn 

/ .OU 




1 nn 


1 fin 

lo.uu 


1 9 nn 


9.95 


8.90 


7 Si^ 

/ .OO 


7 18 


fi 7n 
o. / u 


fi 9n 

O.ZU 




1 nn 


lo.oo 


1 1 1 n 


9.35 


8.35 


7 9n 

/ .ZU 


o.oo 


fi 1 n 

O. lU 


o.oo 


Incr 717 
log Zj/^q 




U.Z04 


n 91^4 


0.254 


0.254 


n 9'^4 

U.ZO^ 


n 91^4 
u.zo^ 


n 9'^4 

U.ZO^ 


n 9'^4 
u.zo^ 








n 91^4 


0.254 


0.254 


n 9'^4 


n 91^4 


n 9'^4 


u.zoo 




n 91^4 


U.Z04 


n 91^4 


0.254 


0.254 


u.zoo 


n 91^"? 
u.zoo 


u.zoo 


u.zoo 




n 91^4 


U.Z04 


n 91^4 


0.254 


0.253 


u.zoo 


u.zoo 


u.zoo 


u.zoo 




n nfiQ 


n nfiQ 
-u.uoy 


n nfiQ 
-U.uoy 


-0.070 


-0.070 


n n7n 
-u.u ( u 


n n7n 

-u.u 1 u 


n n7n 
-u.u / u 


n n7n 

-u.u 1 u 




n nKQ 


n nfiQ 
-u.uoy 


n n7n 
-u.u/u 


-0.070 


-0.070 


n n7n 
-u.u ( u 


n n7n 
-u.u/u 


n n7n 
-u.u ( u 


n n7n 
-u.u/u 




n nfio 
-u.uoy 


-u.uoy 


n n7n 

-U.Ui u 


-0.070 


-0.070 


n n7n 
-u.u ( u 


n n7n 

-U.Ui u 


n n7n 
-u.u ( u 


n n7n 
-u.u ( u 




-u.uuy 


-u.u / u 


-u.u ( u 


-0.070 


-0.070 


-u.u / u 


n n7n 
-u.u ( u 


n n7n 

-u.u ( u 


n n7n 
-u.u ( u 


[Mo- /Pol 


n 4fi7 


U.'iU / 


n 4fi7 


0.467 


0.467 


n 4K7 


n 4fi7 

U.^O ( 


n 4K7 


n 4R7 

u.^o / 






U.4D / 


n AP.7 


0.467 


0.467 


n /IK7 


n AR7 


n /IK7 


n /1K7 




n 4fi7 


Af{7 
U.4U / 


n 4fi7 


0.467 


0.467 


n 4R7 


n 4fi7 


n 4R7 


n 4R7 

U.^0 / 




n 4fi7 


/Lfi7 


n 4fi7 


0.467 


0.467 


n 4R7 


n 4fi7 


n 4R7 


n 4R7 


C — JL 


9 1 (?7 


Z.104 


9 1 4n 


2.126 


2.114 


9 nQQ 

z.uyy 


9 ns9 

Z.UoZ 


9 nfis 

Z.Uuo 


9 n^i/i 

Z.U04 




9 1 fi7 


9 147 


9 1 9Q 


2.112 


2.089 


9 nfis 

Z.UUo 


9 n47 
Z.U^ 1 


9 n'^1 

Z.UOl 


9 ni 9 




2.167 


2.139 


2.110 


2.074 


2.043 


2.013 


1.988 


1.963 


1.938 




2.167 


2.133 


2.098 


2.058 


2.023 


1.989 


1.961 


1.932 


1.904 




1.633 


1.627 


1.621 


1.614 


1.609 


1.603 


1.596 


1.590 


1.584 




1.633 


1.624 


1.615 


1.608 


1.599 


1.590 


1.581 


1.574 


1.566 




1.633 


1.620 


1.607 


1.593 


1.579 


1.567 


1.553 


1.539 


1.525 




1.633 


1.618 


1.603 


1.586 


1.570 


1.553 


1.538 


1.522 


1.506 


^06 — -^814 


0.577 


0.575 


0.573 


0.570 


0.569 


0.566 


0.564 


0.561 


0.560 




0.577 


0.574 


0.571 


0.569 


0.565 


0.561 


0.559 


0.558 


0.563 
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Table 2 — Continued 



logr/re -1 -0.875 -0.75 -0.625 -0.5 -0.375 -0.25 -0.125 

0.577 0.573 0.568 0.562 0.559 0.562 0.557 0.547 0.537 
0.577 0.572 0.566 0.560 0.559 0.557 0.546 0.534 0.524 

Note. — The first, second, third and forth rows of each physical quantity correspond to those for 
the color gradient of A{B — i?)/Alogr = —0.05, —0.07, —0.11, and —0.13 mag/dex, respectively. 
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Table 3: Sample Elliptical Galaxies 



ID^ 








(^606 — hl4)AB 




B/T-" A(V606 - /814)ab'' 






(mag) 


(mag) 


(mag) 


(kpc) 




(mag) 


2-456.0^ 


0.089 


18.28 


-19.20^ 


0.610 


1.4 


- 


-0.02 


4-241.1 


0.318 


20.33 


-20.94 


0.422 


3.1 


0.66 


- 


2-121.0 


0.475 


20.03 


-20.94 


1.071 


3.5 


0.87 


0.10 


3-790.0 


0.562 


20.93 


-21.92 


1.152 


2.8 


0.63 


0.21 


3-321.0 


0.677 


20.92 


-22.17 


1.490 


7.0 


0.81 


0.09 


4-744.0 


0.764 


20.41 


-23.12 


1.609 


5.1 


0.79 


0.06 


4-493.0 


0.847 


21.17 


-22.89 


1.670 


4.0 


0.81 


0.02 


2-251.0 


0.960 


20.86 


-23.65 


1.055 


7.9 


1.00 




4-752.0 


1.013 


20.88 


-24.10 


1.855 


6.9 


0.87 


-0.14 


4-928.0 


1.015 


21.75 


-23.03 


1.347 


3.7 


0.81 


-0.4 



'^Williams et al. (1996). 
^Cohen et al. (1996). 

■^Abraham et al. (1999). We converted the Mp^^qw^ab magnitudes into those in the adopted 
cosmology. 

'^Bulge-to-total luminosity ratio in the /8i4-band taken from Marleau & Simard (1998). 

•^Shifted zero-points to compare the model color gradients for typical elliptical galaxies with 
observed ones (see §4). 

•''The B/T ratio of this object is not in Marleau & Simard (1998). However, we add it to 
our sample because its spectrum (from the Hawaii active catalog; Cowie 1997) as well as 
morphology show this galaxy is clearly a typical E or SO. 

^Bouwens, Broadhurst, & Silk (1998). 

'^fej-band absolute magnitude calculated by Bouwens et al. (1998). We also converted the 
magnitude into that in the adopted cosmology. 
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Figure Caption 

Fig. 1 — Veoe — hu colors of our sample galaxies in a 10 kpc aperture as a function of 
redshift. Solid line shows a color track for a passively evolving elliptical and dotted line for a 
non-evolving elliptical. Filled circles indicate red galaxies whose colors are nearly consistent 
with that of the passively evolving or non-evolving ellipticals. Open circles indicate the 
remaining bluer galaxies. 

Fig. 2 — Azimuthally averaged radial surface brightness profiles of our sample galaxies 
in the /8i4-band. (a) Red galaxies of which colors are consistent with that of a passively 
evolving galaxy, (b) Galaxies of which colors are bluer than that of a passively evolving 
galaxy (see also Figure 1). Filled circles show the data points and solid lines show the fitted 
lines with the r^/^ law. 

Fig. 3 — Veoe — -^8i4 color maps for (a) the red galaxies whose colors are consistent with 
that of a passively evolving galaxy, and (b) the "blue" eUipticals. Plus sign of each object 
denotes the centroid of the galaxy in the /8i4-band image. The circle represents an effective 
radius centered on the centroid. Object ID, redshift, and effective radius are also shown. 
Note that a color bar is different from galaxy to galaxy, to show the color distribution 
clearly. 

Fig. 4 — Observed color gradients of the red elliptical galaxies together with the model 
gradients by the age gradient. Filled circles represent the data points. Crosses refer to 
the data points which are not used in deriving slopes of color profiles. The models show 
the predicted color gradients seen at each object's redshift. Solid lines correspond to 
A{B — i?)/Alog(r/re) — —0.09 ± 0.02 mag/dex at ^ = 0, and the dotted lines correspond 
to A{B - i?)/Alog(r/re) = -0.09 ± 0.04 mag/dex at z = 0. Zero-points of the observed 
colors are slightly shifted by A(V606 — -^814) in Table 3 to compare the observed gradients 
with the models. 
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Fig. 5 — Same as Fig. 4, but for the case of the metaUicity gradient. 

Fig. 6 — Vqo6 — Isu color gradient versus redshift for our red sample galaxies. Filled 
circles represent color gradients of the sample ellipticals, which are determined within 
effective radii without the innermost regions of the galaxies. Sohd and dotted hues indicate 
the evolutions of color gradients caused by the metaUicity gradient and the age gradient, 
respectively. (See text for the model.) 

Fig. 7 — Color gradients of the galaxies which have relatively blue colors. The color 
gradient of 4-928.0, which has a rather steep color gradient, is plotted together with a 
model color gradient made by age gradient fitted to the observed one by tuning ages of 
stellar populations in the model galaxy. (See text for the model.) The observed color for 
4-928.0 is shifted by 0.4 mag to compare the model gradients. 




ure 2(a) 




ure 2(b) 




2-121.0 (z = 0.475, rg=3.5 kpc) 4-493.0 (z = 0.847, rg=4.0 kpc) 




X (pixel) X (pixel) 




3-321.0 (z = 0.677, r =7.0 kpc) 




* 1 pixel = 0.04 arcsec 



20 40 60 80 100 



X (pixel) 



ure 3(a) 



4-241.1 (z = 0.318, r =3.1 kpc) 4-928.0 (z=1.015, r =3.7 kpc) 




X (pixel) X (pixel) 



2-251.0 (z = 0.960, t^ = 7.9 kpc) 




1 pixel = 0.04 arcsec 



20 40 60 80 100 



X (pixel) 



ure 3(b) 



Age Sequence 



0.7 



I 0.6 

to 
o 

CD 

> 



0.5 



0.7 



I 0.6 



CD 
O 
CO 



> 



0.5 



1.3 



I 1.2 



o 

CD 



1.1 



1.5 



I 1.4 



CD 
O 



> 



1.3 



"T 1 1 1 1 1 1 1 1 1 1 1 r 

z= 



H — — \ — \ — \ — \ — i — \ — \ — \ — \ — \ — ^ 

-1 -0.5 
z= 0.089 2-456.0 




H \ \ \ \ \ \ \ \ \ h 

-1 -0.5 

z= 0.475 2-121.0 




H — — \ — \ — \ — I ■ n I — \ — — \ — I- 

-1 -0.5 
z= 0.562 3-790.0 




1.7 



I 1.6 

CO 

o 

CD 



1.5 



1.8 



I 1.7 

CD 
O 
CO 

> 

1.6 



1.8 



I 1.7 

CO 

o 

CD 



1.6 



1.8 



CD 
O 



> 




I 1.7 



1.6 



z= 0.764 4-744.0 




-1 -0.5 
z= 0.847 4-493.0 




■1 -0.5 
z= 1.013 4-752.0 








-0.5 



-0.5 



log r/r^ 



log r/r^ 



Metallicity Sequence 



0.7 



I 0.6 

to 
o 

CD 

> 



0.5 



0.7 



I 0.6 



CD 
O 
CO 



> 



0.5 



1.3 



I 1.2 



o 

CD 



1.1 



1.5 



I 1.4 



CD 
O 



> 



1.3 



"T 1 1 1 1 1 1 1 1 1 1 1 r 

z= 



H — — \ — \ — \ — \ — i — \ — \ — \ — \ — \ — ^ 

-1 -0.5 
z= 0.089 2-456.0 




H — — \ — \ — \ — \ — — \ — \ — \ — \ — \ — ^ 

-1 -0.5 
z= 0.475 2-121.0 




H \ \ \ \ \ \ \ \ \ 1- 

-1 -0.5 

z= 0.562 3-790.0 




1.7 



I 1.6 

CO 



CD 
O 
CO 



> 



o 

CD 



> 



1.5 



1.8 




I 1.7 



1.6 



1.8 



I 1.7 



1.6 



1.8 



I 1.7 

CD 



1.6 



z= 0.764 4-744.0 




H \ \ \ \ \ \ \ \ \ h 

-1 -0.5 

z= 0.847 4-493.0 




■1 -0.5 
z= 1.013 4-752.0 




-0.5 



-0.5 



log r/r^ 



log r/r^ 




log r/r^ 



